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CALCULATION OF ELECTRON CONCENTRATION 
FOR A BLUNT BODY AT ORBITAL SPEEDS AND COMPARISON 


WITH EXPERIMENTAL DATA 

By Charles J. Schexnayder, Jr., Paul W. Huber, 
and John S. Evans 
Langley Research Center 

SUMMARY 

Calculation of electron concentration for a blunt-nosed body at orbital speeds has 
been made by use of pure-air nonequilibrium chemistry throughout the flow field and for 
the same flow field with alkali contamination (in local ionic equilibrium) in the boundary 
layer. Plasma diagnostic data from onboard Langmuir probes, reflectometers, and 
antenna voltage standing wave ratio (VSWR) , along with passive data from attenuation of 
beacon and telemetry signals, are available for the first two RAM C reentries. These 
data are compared with calculated flow-field ionization and radio -frequency (RF) signal 
attenuation for both inviscid flow and for inviscid flow merged with boundary layer. Good 
agreement between plasma theory and data is obtained in the altitude range between about 
30 km and 71 km. At altitudes greater than 71 km, ambipolar electron-ion diffusion and 
wall recombination are shown to have a large effect on peak electron concentration in the 
plasma. For RAM C-I at altitudes less than 60 km, there is a large increase in electron 
concentration in the boundary layer due to ionization of alkali impurities from the heat 
shield. Indications of boundary-layer transition along the body and of ablation rate on 
the nose were obtained from the low-altitude reflectometer data. 

INTRODUCTION 

The free electrons which are produced in hypersonic-flow-field plasmas are the 
dominant factor in the entry radio blackout phenomenon. Therefore, in order to define, 
circumvent, or alleviate the blackout problem, it is important that one be able to deter- 
mine the magnitude and distribution of electrons in the plasmas near communication 
antennas, as well as to determine the interactions which occur between the plasma, the 
antennas, and the propagating electromagnetic signals. Theoretical evaluation of the 
electron density levels and profiles for a given entry communications problem requires 
precise and detailed consideration of complex systems of chemical and thermodynamic 
processes for the complete shock layer, since free electrons are present only as trace 



species in time -dependent systems. (See ref. 1.) Even with careful analysis, uncer- 
tainties in the flow species to be considered, the reaction rate coefficients, viscous inter- 
actions, antenna near -field effects, and so forth combine to make applicability to many 
real-life situations in doubt, particularly those involving ablation contamination (ref. 2) 
or high-altitude blackout. 

The problem of radio blackout has led to the establishment of a flight project at the 
Langley Research Center called the Radio Attenuation Measurements Project (RAM). The 
basic objectives of Project RAM have been to obtain a better understanding of the factors 
which influence transmission of radio waves through entry vehicle plasmas and to search 
for ways to reduce or eliminate blackout. 

The RAM project has included three series of flight payloads; namely, RAM A, 

RAM B, and RAM C. All were spherically blunted 9° half-angle cones with nose diameters 
of 5.08 cm (2 inches), 20.32 cm (8 inches), and 30.48 cm (12 inches), respectively. Maxi- 
mum velocity for both the A and B series was about 5.5 km/sec (18 000 ft/sec) whereas 
the maximum speed for the C series was about 7.6 km/sec (25 000 ft/sec). Descriptions 
of the flights and the results were published in references 3 to 8 and references 9 to 12 
dealt with an analysis of flight results. 

The purpose of this paper is to compare flight results from RAM C-I and C-II with 
theoretical calculations and to draw some conclusions from the comparison. In addition 
to the reflectometer data of RAM C-II (ref. 8), this paper uses attenuation and signal 
blackout data from the receiving stations (refs. 7 and 8) to infer the plasma flow-field 
properties. Also, each RAM C flight carried a rake of Langmuir probes (refs. 7 and 13) 
which recorded ion concentration profiles on the aft part of the payload. 

Included in the paper are appendix A which deals with the clean-air flow field, 
appendix B which deals with a mass conservation computer program, appendix C which 
gives details of the boundary -layer treatment, appendix D which shows the effects of 
ablation on the electron concentration in the boundary layer, and appendix E which gives 
the interpretation of low-altitude reflectometer data. 

SYMBOLS 


A 

Cp 

C P 


surface area, cm^ 

p - p 

normalized pressure coefficient, 

p ns " p °° 

specific heat at constant pressure of air mixture, erg-g”^-K 
reaction energy divided by RT 0 for Zth reaction 
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d N 
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Hi 

K l 
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l 

M 

m 

m 

N 

n a 

N e 

NRe 

P 

R 


body nose diameter, cm (in.) 

fraction (by weight) of alkali contaminant in heat shield 

signal frequency, MHz 

degeneracy factor 

enthalpy of ith species, erg-mole"^ 

equilibrium constant for 7th reaction, mole-cm - ^ or dimensionless 
reaction rate constant for 7th reaction, cm -mole" -s or cm -mole" -s 
rarefaction parameter 
reaction identification number 

molecular weight, g-mole"*; also Mach number; also third body in a 
chemical reaction 

flow rate, g-s - ^ 

-2 -1 

ablation gas surface flux, g-cm -s 
quantity used in Saha equation (see eq. (13)) 

Avogadro constant, 6.02217 X 10^ particles -mole 

O 

electron concentration, electrons-cm 
Reynolds number 

o 

pressure, dyne -cm 

2 -5 

contribution of 7th reaction rate to rate of change of ith species, mole -g 
gas constant, erg-mole"^-K"^ or cm^-atm-mole~^-K"^ 
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distance from body axis to body surface, cm 


r b 

r s 

(SL) a 

(SL) t 

S 

s 

T 
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u 

V 

Xi 

x 

y 

z 

a 

r p 

y 
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distance from body axis to shock front, cm 
signal loss due to absorption, dB 
total signal loss, dB 
entropy, erg-mole - * 

distance along streamline from shock front, cm 

temperature, K 

time from lift-off, s 

flow velocity, cm-s - ^ 

effective ionization energy, eV 

mole fraction of alkali species 

distance from nose along body axis, cm 

distance from body along normal to body surface, cm 

compressibility factor 

absorption coefficient or fraction of ionizable species singly ionized 
power reflection coefficient 
specific heat ratio 

specific concentration of electrons, mole-g - * 

specific concentration of ith species, mole-g -1 

value of y at shock front for a given value of x/d^, cm 
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distance to edge of boundary layer, cm 


-'m 


6 * 


effective boundary-layer thickness for merged profile, cm 


boundary-layer displacement thickness, cm 


e angle between axis and tangent to shock wave , deg 

V^Yi, mole-g~l 


6 b 


V- 




CO 


angle between axis and tangent to spherical nose, deg 
angle between axis and conical body surface, deg 
viscosity, g-cm"^-s"^ 

reduced mass of colliding species for Ith reaction, g-mole" 
collision frequency of an electron with other particles, s~^ 
density, g-cm"^ 

cross section for Ith reaction, cm 2 
angular frequency, 27 rf, radian-s”^ 


Subscripts: 


abl 


alk 


ablation material 


alkali material 


cr 


body 


critical value 


cone 


electron or quantity at boundary-layer edge 



f endothermic rate direction 

i species identification number 

ns value at shock front on stagnation streamline 

o standard conditions (p Q = 1 atm, T 0 = 273.2 k) 

p peak 

r exothermic rate direction 

s value at shock front on a general streamline 

°o ambient flight conditions 

1,2, 3, 4 reflectometer body stations (see table I) 

Special notations: 

cr calculated N e = N e>cr 

cr* Langmuir probe measured N e = N e>cr 

tr transition from laminar to turbulent flow in boundary layer 

Primes denote transformed coordinates and circumflex, a normalized quantity. 

RAM C FLIGHT DESCRIPTION 

There were three RAM C flights but only the first two are discussed here, since the 
last one took place after this paper was written. The RAM C payload was a 9° half-angle 
blunted cone with a nose diameter of approximately 30.5 cm and entered the earth's atmo- 
sphere at a velocity of 7.6 km/sec. Figure 1 gives the trajectory of the first two flights. 
RAM C-I (ref. 7) was a material-addition experiment with water injection at the nose and 
at side positions. Four communications antennas were also used for passive diagnostic 
instrumentation. On RAM C-H (ref. 8), the principal effort was on a plasma diagnostic 
experiment based on 15 reflectometers. There were also four communications antennas. 
Each payload carried a rake of eight Langmuir probes which extended 7 cm into the aft 
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flow field in a direction normal to the surface. Table I lists the diagnostic instruments 
and their locations for the two flights. During the entry data period, the payloads had a 
spin rate of about 3 rps and a near -zero angle of attack, except for a few short periods 
where it was as high as 5°. 

It might be noteworthy to mention that although both payloads had the same con- 
figuration, there were two different ways of protecting the payload from the heat generated 
at the nose. (See sketch in fig. 1.) In the RAM C-I flight, NARMCO 4028, a phenolic 
graphite charring ablator, was used at the nose. In the RAM C-H flight for the primary 
data period 100 to 56 km, a beryllium (Be) nose cap positioned over a teflon ablator 
allowed reflectometer data to be taken in a clean-air flow. The nose cap was ejected at 
approximately 56 km; the teflon remained as an ablator at lower altitudes. Both payloads 
used teflon on the afterbody for heat protection. 

An analysis of a sample of the NARMCO 4028 used to fabricate the heat shield 
showed that it contained about 1100 pg/g sodium (Na). The analysis was incapable of 
detecting potassium (K), if present, at less than 3600 p g/g which means that the ablator 
could have contained up to 4700 pg/g of alkali. The Be nose cap and the teflon were found 
to be free of any significant amounts of alkaline impurities. For a sample of teflon the 
analysis showed the alkali content to be less than 5 ppm. 

THEORY 

Clean-Air Flow Field 

The techniques described in reference 1 provide a good approximation to the struc- 
ture of the flow field around a spherically blunted cone such as the RAM payload. Since 
a knowledge of the electron concentration in the flow field is the major factor in any study 
involving radio blackout, the ability to calculate the concentration and distribution of the 
electrons is important. As in reference 1, emphasis is placed on the chemical kinetics 
and the influence of the boundary layer in computing the properties within the blunt-body 
flow field. The reader is referred to appendixes A, B, and C for a detailed description 
of the flow-field calculations. 

In the inviscid flow part of the calculations, density, velocity, temperature, and gas 
composition are computed as functions of distance along a streamline. Velocity and 
density are then used in a mass balance program to obtain y-coordinates for the stream 
tubes at prespecified x/d N locations. At altitudes greater than 30 km, boundary-layer 
corrections to inviscid results are obtained by following streamlines into the laminar 
boundary layer. The stream-tube method used accounts for the effects on the flow field 
of boundary-layer displacement, streamline swallowing, and vorticity interaction normal 
to the flow direction. The effects of finite-rate chemistry on a streamline position and 
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shock shape are also included. A shortcoming of the method (even with boundary-layer 
corrections) is an inability to include the effects of electron-ion diffusion. This effect, 
as will be shown later, is not negligible for N e profile peak values for RAM C at alti- 
tudes greater than 70 km. 


Contaminated Boundary Layer 

Ablation from the RAM C-I heat shield of materials known to contain easily ionized 
impurities such as sodium (Na) and potassium (K) made it necessary to account for an 
expected increase in electron concentration in the boundary layer. Appendix D gives the 
details of the calculation along with typical electron concentration profiles useful for com- 
paring with RAM C-I flight data. Alkali ionization is assumed to be decoupled from the 
air-ablation chemistry and is in equilibrium with the local nonequilibrium thermal proper- 
ties of the air boundary layer as determined previously. No calculations have been made 
for the air-teflon boundary layer, for reasons to be discussed later. 

RESULTS 

Comparison of Data With Plasma Calculations 

Values at the peaks of pure air N e profiles (see appendix A) are plotted in fig- 
ure 2 as a function of altitude and velocity at body locations appropriate for comparison 
with data from the RAM C-II flight. Electron concentrations, as derived herein from 
reflectometer data, are plotted as horizontal bars at values of N e>p = N e cr . (The low- 
altitude bars are dashed to indicate the presence of uncertainties to be discussed later.) 
The length of a bar represents the period during which the reflection coefficient changed 
rapidly. Examples of the region of rapid change in the reflection-coefficient data are 
shown in figure 3, where reflection -coefficient records from reference 8 are shown with 
arrows denoting the times for theoretical N 6jCr , Langmuir-probe-measured N ejCr , and 
where transition from laminar to turbulent flow probably occurred in the boundary layer. 
The region where body motion was observed is shown as well as the region of fast abla- 
tion. (See appendix E and ref. 8.) Figure 4 shows samples of reflectometer phase data 
similarly marked. Shock-tube tests of the X- and S-band reflectometers in reference 14 
indicated that critical N e was associated with the end rather than with the beginning of 
the rise to high reflection coefficient (that is, the low-altitude high r p ends of the bars) 
for the collision frequency range (v/io < 1.0) appropriate to the onset of critical N e . For 
high -collision frequency (v/u 2 1.0), the theoretical indications are that N e cr occurs 
at low values of reflection coefficient; thus, for the low-altitude reflectometer decay data, 
the N e>cr point is expected to be at the low-altitude low r p end of the bars. 
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In addition to the reflectometer data, the rise of VSWR in the VHF antenna feed 
circuits is shown in figure 2. Langmuir probe data for the RAM C-II flight (ref. 13) are 
shown as a cross-hatched area representing envelopes (due to spin and angle of attack) 
of data for the outermost probe (y = 7.08 cm). (The highest ion concentration values mea- 
sured on each profile were always at the outermost probe.) That part of the Langmuir 
probe data which is not hatched indicates the altitude range for which probe temperatures 
were above 812 K. These data may have been affected by thermal degradation of the 
beryllium oxide insulation. (See ref. 7.) The VHF pattern data in figure 2 represents 
the altitudes at which the antenna patterns (seen because of vehicle spin) were observed 
briefly to be circular. (See ref. 15.) The analysis in reference 16 for a similar antenna 
array indicates that the occurrence of a circular pattern is a definite and sharp indicator 
for the condition N e ,p = N e cr . This observation is an excellent check on the other 
diagnostic data. Signal attenuation data for X-band, C-band, and VHF frequencies for 
the RAM C-H flights are converted to N e? p values for display in figure 2 by using the 
plane-wave nonhomogeneous plasma model given in reference 17. This conversion is 
done by first assuming a generalized profile shape (as a function of altitude) and arbi- 
trarily letting the N e magnitude vary. Values of N e> p are then found whose profiles 
yield theoretical attenuations equal to the measured attenuations. The N e profile shapes 
assumed in this method are similar to those found in appendix A. 

Values at the peaks of the air plus alkali N e profile curves (see appendix D) are 
plotted in figure 5 as a function of altitude and velocity for x /d^j of 2.0 and 4.0 and are 
compared with data from the RAM C-I flight. Curves representing pure air N e> p are 
also shown. Langmuir probe data (ref. 7) are shown as a cross-hatched area representing 
envelopes of data for the outermost probe. In addition to N e p values derived from 
attenuation data for X-band, C-band, and VHF frequencies, the rise of VSWR in the VHF 
antenna feed circuits along with the VHF pattern data are also shown. 

Inspection of the theory and flight results in figures 2 and 5 for the altitude regime 
greater than 60 km brings up some interesting points. (1) There is generally good agree- 
ment among the various diagnostic techniques for this altitude range. (2) In view of the 
agreement of RAM C-I and C-H data, there does not appear to be any effect of alkali 
ionization at altitudes greater than 65 km. (3) At all body stations and in the altitude 
range above 71 km, calculated peak N e values are seen to be much higher than the mea- 
sured values. However, the theoretical model utilizes the stream-tube concept and hence 
does not allow for inclusion of effects of electron-ion diffusion or for wall recombination 
effects. 

Estimated correction factors for these effects are obtained by correlating RAM C 
data with shock-layer solutions from the literature which include diffusion and wall effects. 
The results are plotted in figure 6. The hatched region is one which includes theoretically 
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and experimentally determined N e> p reduction factors taken from results of refer- 
ences 18 to 22. This factor is the one by which N e>p for an inviscid profile is reduced 
because of the effects of viscosity, ion diffusion, and wall recombination. In order to plot 
the flight data in this figure, they are divided by N e p values from inviscid curves. All 
the reference data are for the nose region of blunt-nosed bodies with the exception of 
reference 19 which is for x/d N = 4. Values appropriate to the stagnation point are gen- 
erally near the high-altitude side of the hatched region, and values for the nose -body 
juncture are near the low-altitude side. The reference data are correlated on the RAM C 
altitude scale on the basis of the rarefaction parameter (used to scale results in ref. 21) 


k 2 - L& ' 1 ^°° dN 

2ye© 2ji ns 


( 1 ) 


The correlation is seen to be indicative of the effect of ion diffusion and wall effects on 
the RAM C shock-layer peak N e . The effects begin at about the altitude predicted 
(k2 « 100). 

The comparison between air theory and RAM C-n flight data in figure 2 is reason- 
ably good at altitudes greater than 30 km but less than 60 km. The low-altitude ends of 
the bars for the X3 and X4 reflectometer data are believed to be near the critical electron 
concentration. (See ref. 14.) 

A comparison between the RAM C-I flight data (30 km < altitude < 60 km) and the 
two sets of theoretical curves (representing air and air plus alkali theories) in figure 5 
indicates that the data from the RAM C-I flight correlate better with the contaminated 
boundary-layer model than with the air model. This condition is true for the lower alti- 
tude part of the flight data. Above 50 km, the data drop off where the air plus alkali theory 
is rising. This drop-off of N e> p is probably due to failure of the local equilibrium 
assumption in this region (probable delay in production of alkali ions). This figure illu- 
strates that there is a large increase in the electron concentration in the boundary layer 
due to ionization of alkali impurities from the heat shield of RAM C-I for altitudes less 
than 60 km. Comparison of flight data (altitude greater than 30 km) in figures 2 and 5 
indicates that alkali ionization is probably responsible for the difference in the electron 
concentration between the RAM C-I and RAM C-II flights. 

For comparison of the low-altitude flight data and theory (altitude less than 33.7 km), 
a velocity scale is used in figures 2 and 5. Table III gives the velocity-altitude pairs 
used. For the RAM C-I flight, N e>p values were obtained from attenuation data for 
X-band, C-band, and VHF frequencies by a method previously described. For the 
RAM C-II flight, N 6j p values were obtained from attenuation data for VHF frequencies 
only, since signal recovery data for X- and C-band were not obtained because of tracking 
problems. Reflectometer data from RAM C-II are represented in figure 2 by bars as 
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previously described. The theoretical curves shown in figures 2 and 5 are computed 
for an air model where any electrophilic effects of teflon ablation into a laminar boundary 
layer are not considered since the boundary layer is thin relative to the plasma, and 
reduction of electrons would not appreciably affect N e ,p. However, calculations are 
made for a model containing the effects of alkali ablation in a laminar boundary layer 
(appendix D) since the shock layer N e?p is thus greatly increased. 

The data from the RAM C-I flight in figure 5 agree reasonably well with the theo- 
retical curves labeled "air plus alkali," which indicates that the N e>p values are two to 
three decades higher than the N e>p values computed by use of an air model. This dif- 
ference is attributed to alkali impurity ionization in the boundary layer. The Ne,p 
values obtained from the VHF attenuation data for the RAM C-II flight agree fairly well 
with the theoretical air curves in figures 2(d) and 2(e). The reflectometer data bars 
(fig. 2) do not agree with the theoretical curves. They show a much steeper drop-off of 
electron concentration than does theory as the payload slows down. 

Uncertainties in Theory and Data 

Since the low-altitude reflectometer data do not agree well with the plasma theory 
or deduced VHF data, possible uncertainties in both theory and data should be examined. 
One source of uncertainty in the theory is in the values of reaction rates used to deter- 
mine plasma composition. Theoretical curves for pure air are given in figure 7 along 
with the effects of changing or deleting certain reaction rates. Changing the 15 dissocia- 
tion reaction rate constants (1 to 15, table H) by factors of 10* and lO - * produces N e>p 
changes as indicated by the hatched region for three of the body locations. The factors 
for the rate changes are arbitrarily taken and are not intended to indicate that the dis- 
sociation rates of 02> N 2 , and NO are uncertain by such large factors. Resultant changes 
in the N e ,p level are not large. This result would be expected since reference 23 
reported similar equilibrium and nonequilibrium N e values for the nose flow at 24.5 km 
(Velocity = 5.12 km/sec) and indicated that N e was not sensitive to the dissociation rate 
of O 2 , N 2 , and NO at this altitude. Deleting all reactions which involve O” and Og” causes 
the electron density to be larger at velocities below 4.5 km/sec for an aft station. Fac- 
tors of 10* and 10“1 for the rate constants involving O" + M rz + e" + M and 
Og - + M ~ Og + e" + M gave a change in the N e>p of only 20 percent at a velocity of 
3.9 km/sec for the aft body stations. At higher velocities the changes in the rate constant 
for O" and O 2 ” do not produce any noticeable changes in the N e p values. The exclu- 
sion of NO 2 - and reactions involving this molecule may cause some error in the low- 
speed calculations. Reference 24, which has equilibrium results for air, reports that 
NO 2 ” becomes a dominating factor in determining the N e at pressures higher than 
1 atmosphere and temperatures less than 4000 K. However, the combination of pressure 
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and temperature where the reflectometer data are obtained does not match the conditions 
needed to make NOg" important. 

One way the air calculation could be appreciably perturbed would be by the combi- 
nation of large boundary-layer thickening, mixing, and chemical action of the teflon abla- 
tion products. The thickening and mixing would presumably occur by the combined effects 
of transition, increased mass injection, body motion, surface roughness, and so forth. 
These effects would have to thicken the boundary layer on the nose to «10 times the com- 
puted laminar value for ablation products to reach the required point (outermost point 
where N e ~ N 6jCr ) in the inviscid plasma for agreement with reflectometer data. Such 
large thickening seems unlikely since turbulent boundary- layer calculations show that 
transition would only cause about a factor of 2 increase in thickness. The maximum 
blowing parameter for the teflon is theoretically ^0.06, but the effects of blowing have not 
been determined. Also, material injection studies have shown it to be difficult to cause 
even liquids under pressure to penetrate such high density flows; therefore, gaseous abla- 
tion products would not be expected to penetrate the flow so deeply. In addition, any large 
thickening of the boundary layer would move the inviscid plasma out even further because 
of displacement effects. 

The ablation of teflon (ref. 25), a polymer of tetrafluoroethylene (C2F4), occurs, 
sublimation producing the monomer gas. The high -temperature behavior of teflon is 
discussed in reference 26. Decomposition and oxidation of C2F4 in the hot air flow will 
produce a wide variety of products, some of which are CO, CO2, F, CF4, CF 2 0, CF3, CF2, 
and CF. (See refs. 27 to 29 for the chemistry involved.) The structure of the laminar 
ablating air -teflon boundary layer is given in reference 30. 

If the teflon-air products could somehow deeply penetrate the nose-plasma layer, 
the chemical action required to fit the observations must be either electron attachment 
or interference with electron production. The computed ablation mass is theoretically 
about 6 percent of the mass within the required plasma stream tubes. Reference 31 has 
computed only about 1 percent of the electrons attached for a 5-percent teflon-air equi- 
librium mixture at about the temperature (6000 K) and density of RAM C-n nose flow in 
this low-altitude flight regime. Reference 32 has reported no difference in the conduc- 
tivity between air and air -teflon mixtures (“3 percent teflon in air) behind normal -shock 
waves where the gas temperature was 3200 K to 6500 K and density on the order of 0.1p o 
to 0.2 p Q . It has been shown experimentally in reference 33 that for mixtures of 
C^F^-air above 3000 K, electron quenching was not obtained. 

The agreement of RAM C-I low- altitude data with the alkali theory tends to rule out 
any large teflon electrophilic action in the boundary layer (x/d^- > 0.8) since the body 
coating injects more teflon into the boundary layer than there are NARMCO products 
present due to nose ablation. With the teflon present and if attachment of electrons 
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occurred, then the electron concentration for RAM C-I would not show a large spike 
(appendix D) in the boundary layer and attenuation data would not agree with the air plus 
alkali theory. 

Although theoretical calculation of attenuation for the RAM C VHF case is not avail- 
able, some feeling for the possible uncertainties in the interpretation of the VHF attenua- 
tion data can be had by examining figure 8 where theoretical plane wave VHF attenuation 
curves for 226 MHz at x/d^j = 4 are compared with recovery data at x/d^ = 3.8 for 
the RAM C-H flight. A curve showing attenuation due to absorption only (a lower limit 
for plane-wave signal loss) is shown along with a theoretically computed total signal loss 
curve by using the plane-wave nonhomogeneous plasma model in reference 17. The dif- 
ference in the level of attenuation between the two curves is due to plasma reflection 
losses, and to each curve must be added any mismatch losses which may actually be pres- 
ent. (The magnitude of the RAM C VHF mismatch loss is not known.) During the body- 
motion period, the angle of attack of the vehicle changed by as much as 5° so that calcu- 
lations of (SL)^ were made at 40.0 km and 30.6 km by using the N e profiles given in 
appendix E. The hatched region gives the limits of attenuation in the body-motion period. 

Several aspects of figure 8 suggest that the total -signal -loss calculation gives 
results that are preferred over the absorption-only case. (This statement is not to imply 
that plane-wave theory is valid for VHF, but merely that the total RAM-C losses due to 
absorption, reflection, mismatch, etc., happen to amount to something on the order of 
those calculated by using the plane-wave concept.) First, of course, is the good agree- 
ment of the attenuation data with the total-signal-loss values. Second is the fact that 
VHF recovery did not occur much earlier than observed, as should have been the case if 
plane-wave absorption were the only losses. In particular, there should have been a sub- 
stantial recovery lasting for around 4 seconds during the body-motion period. Recovery 
would not be expected in the case of the higher losses since the reduction during the body- 
motion period (although not plotted) would be only to about the 40 dB (margin) level. 

If one assumes that this total plane wave loss is at least a fair approximation, then 
the fact that VHF recovery did not occur earlier would rule against the low-altitude N e p 
indications given in reference 8, where the faired N e p curves for the aft stations (as a 
function of altitude) showed decreasing N 6j p below about 37 km. Such a decrease in 
N e p, along with the increasing v should result in a decreasing attenuation in the region 
below 37 km (where the theoretical is increasing), recovery occurring at least 3 seconds 
earlier than was observed (by 33 km). 

For the VHF recovery data to be compatible with the low-altitude N e p fairing in 
reference 8 for the RAM C-II flight, large losses, over and above the total-plane-wave 
losses, would have to be present because of antenna mismatch, and so forth. This condi- 
tion can be seen from the difference between observed signal loss at recovery given in 


figure 8 and that computed for the reflectometer indicated peak N e profile (by fairing 
the data in ref. 8, and then using plane-wave theory and maximum plasma thickness to 
obtain an upper limit case). 

Some of the factors which lead to uncertainties in the interpretation of the low- 
altitude reflectometer data were pointed out in reference 8 and are also discussed in 
appendix E. The effect of high electron-collision frequency is to make more difficult the 
interpretation of data from reflectometers for which u/oo > 1. Therefore, identification 
of the point, N e> p = N e>cr , from the Tp record is subject to a large uncertainty in some 
cases. Another uncertainty arises because of the long-duration heat pulse of reentry 
(>20 sec) which causes erosion of the reflectometer dielectric covering. For example, 
the record for K-band at station 1 in figure 3 shows large oscillations which correspond 
closely to dielectric thickness changes rather than changes in N e> p. Ablation may also 
be responsible, in part, for the nonideal character of the decay records of S-band at sta- 
tion 1 and X-band at stations 1 and 2. Finally, the sudden decays of reflected power for 
the X-band records as shown in figure 3 are believed to be associated with transition from 
laminar to turbulent flow in the boundary layer. Hence, there is an uncertainty in inter- 
pretation of the X-band data due to possible effects of changes in profile shape and stand- 
off distance. An interesting aerodynamic result from this effect is that the transition 
apparently occurs at a fairly low Reynolds number (Nr 6 ~ 350 000; see appendix E). This 
result may also suggest the reflectometer as an instrument for transition measurement. 

In summary, much of the low-altitude reflectometer decay data plotted in figure 2 
falls well below the plasma theory and the VHF attenuation data. Reasons for the dis- 
crepancy apparently cannot be found in terms of the chemical kinetics, boundary-layer 
thickness, boundary-layer blowing or electrophilic action effects which were used (or 
neglected) in the development of the plasma theory. The VHF attenuation data, although 
subject to many questions, seem to corroborate the plasma theory. Several of the reflec- 
tometer decays are questionable for indication of N e cr because of deleterious effects 
of high collision frequency, ablation, body motion, and possible sensitivity to N e profile 
shape. 


CONCLUDING REMARKS 

For RAM C-I at altitudes below 60 km, the largest values of electron concentra- 
tion appear to have been present in the boundary layer and were the result of alkali (Na 
and/or K) present in the NARMCO 4028 phenolic graphite charring ablator used for heat 
protection on the nose. There was no evidence at any altitude for alkali contamination 
of the boundary layer for RAM C-H, which had a beryllium nose cap (discarded at 
56 km) and was otherwise protected by teflon. 
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Except for the alkali contamination effects mentioned, electron concentrations com- 
puted for clean air gave satisfactory agreement with observations for altitudes between 
36 km and 71 km. Divergence between theory and data above 71 km for both RAM C-I 
and RAM C-II was found to be due to ambipolar electron-ion diffusion and wall recombina- 
tion. Apparent discrepancies between reflectometer results and calculated values of 
peak electron concentration below an altitude of 36 km are attributed, at least in part, to 
effects of high collision frequency, to changes in the thickness of teflon with time, and to 
effects of laminar -to -turbulent boundary-layer transition. Transition is believed to have 
occurred at a local Reynolds number of about 350 000. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 7, 1971. 



APPENDIX A 


CLEAN-AIR FLOW FIELD 

In this appendix techniques for computing ionization profiles for blunt-nosed reentry 
bodies are presented. Emphasis is on careful consideration of the chemical kinetics and 
of the effects of boundary-layer vorticity interaction. The computations are limited to 
orbital entry velocity or less and to the case of negligible ionizable ablation impurities. 
The body considered herein is a 9° half-angle cone blunted with a 30.5-cm-diameter 
spherical nose, and at zero angle of incidence to the entry flight path. 

The method, which includes several machine programs and one hand operation, is 
tailored for the computation of detailed ionization profiles throughout the supersonic part 
of the shock layer at many points along the entry trajectory. Basically, the uncoupled 
stream-tube approach is used, but coupling is achieved, in effect, by iteration of the 
stream-tube and shock-layer programs. A complete reaction system is employed along 
the streamlines both in the inviscid-flow region and in the merged boundary-layer 
inviscid-flow region, but in the latter region the streamline locations and boundary con- 
ditions are changed according to the dictates of the boundary-layer solutions and the hand- 
merging procedure. Although the method used is not fundamentally different from previ- 
ous methods for blunt bodies with regard to use of the stream -tube technique (ref. 34), 
boundary-layer swallowing (ref. 19), and boundary-layer merging (ref. 23), some of the 
innovations employed lead to more precise specification of the ionization profiles. The 
unique features include the use of finite-rate nonequilibrium shock-layer solutions for 
determination of shock shape, for location of inviscid streamlines, and for specification 
of edge conditions for boundary-layer swallowing and merging. 

Definition of Inviscid-Flow Field 

Figure 9 shows the sphere-cone body and a typical streamline path in the sur- 
rounding flow field. A number of methods for finding the details of such a flow field have 
been published in the literature. In the subsonic region either the direct (ref. 35) or the 
inverse (ref. 36) approach is used, and the solution is continued in the supersonic region 
by the method of characteristics. In these methods the gas properties are usually 
assumed to be either in complete thermochemical equilibrium or to have a frozen com- 
position. The computer time and professional skill required to obtain a complete flow- 
field solution at one flight speed and altitude is substantial. When finite-rate reactions 
are included (ref. 37), the required effort goes up tremendously. 

The stream-tube approach to finite-rate flow property determination avoids the 
complexity of the coupled approach by utilizing the streamline paths and boundary 
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conditions along the streamlines as obtained from equilibrium or frozen flow-field solu- 
tions. The error involved in the resulting ionization profiles is usually small, except in 
the nose regions. In the stream-tube approach used herein, streamline positions are 
adjusted for the effects of finite-rate chemistry through a shock-layer mass-balance 
program discussed in appendix B. 


Furthermore, the necessity to calculate an equilibrium or frozen flow field for each 
trajectory point needed can also be avoided, with small error, by utilizing normalized 
plots of existing solutions. The normalizations which were used to construct inviscid 
streamline positions and pressure distributions are shown in figures 10 and 11. They 
were derived from a number of equilibrium and frozen flow-field solutions. (See refs. 36 
and 38.) In the forms shown, these parameters are relatively independent of Mj, 6 C , 
djq’, and gas properties in the hypersonic regime 15 > Mj > 30, in the range of cone 
angles 6 < 6 C < 15, and for the altitude range 20 < altitude < 70 km.* Even though com- 
putations were made only for a 9° half-angle blunt- cone -shaped body, 6 C is larger than 
9° at the higher altitudes where actual cone angles are replaced by larger effective cone 
angles in order to account for the displacement thickness of the boundary layer. 


For a given set of flight conditions pressure distributions were constructed as 
follows: (1) The shock shape (e - plotted against x/d^ and e/e c plotted against 

n v i ' / * /_ 


0 C 2 and the body pressure distribution ^Cp plotted against x/d$ and p/p ( 


plotted 


against 6 C are found, where 8 C and p c are sharp-cone values based on Mj, 

6>c, and an assumed specific heat ratio; (2) At given x/djvj the variation of pressure and 
entropy across the shock layer is found from figure 11. The entropy parameter is used 
only for determining streamline locations in the shock layer by means of an entropy 
identification, which is chosen to be the entropy behind the equilibrium shock at its entry 
point into the shock layer;** (3) The resulting pressure distribution is curve fitted to pro- 
vide pressure variation along streamlines (fig. 11). 


A standard network of flow-field points composed of intersections between 20 num- 
bered streamlines and 20 body station perpendiculars was devised to provide convenient 
cross-reference points between computer programs and as an aid in plotting results. At 
the beginning of work for given flight conditions, an initial shock line is plotted. (This 
line is subject to later modification by mass-flow considerations.) Each of the numbered 
streamlines originates at the intersection of a body station line with the shock line. Ini- 
tial conditions on a streamline are obtained by iterating the Rankine-Hugoniot shock 


*The upper altitude limit varies with body scale. The limit shown is appropriate 
for djq = 30.5 cm. 

**This entropy value is not used for any other purpose because of the nonequilibrium 
nature of the flow. 
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equations for a gas of free-stream composition and fully excited internal energy across 
an oblique shock tangent to the shock envelope at the entry point. 


Integration Along Stream Tubes 

The next step is to solve the quasi-one -dimensional fluid-flow equations and the 
finite-rate chemical equations by means of step-by-step integration. These equations 
are written as follows: 


pRT/1 dg\ ^s/1 dT\ V V> RT ( 1 dy j \ = 1 dp 
u s 2 \ p s ds / P s u s 2 \ T s ds / P s Us2 ' ?s ds ' P s u s 2 ds 



where 6y is the Kronecker delta. This way of formulating the problem is essentially 
that described in reference 39, in which the solution of the equations is also discussed. 
The first and second equations are derivatives of the equation of state and the energy 
conservation equation, respectively. Substitution of the momentum conservation equa- 
tion — = eliminated the use of the flow speed. The remaining equations give the 

ds pu ds 

production rates of the 13 species. 

The composition of the ambient air is assumed to be O 2 (21 percent) and N 2 (79 per- 
cent) . The species that are considered in the flow are 


Species 

Neutral 

Species 

Positive 

Species 

Negative 

no. 

species 

no. 

species 

no. 

species 

1 

°2 

6 

NO + 

11 

e~ 

2 

n 2 

7 

°2 + 

12 

°2 _ 

3 

O 

8 

N 2 + 

13 

O' 

4 

N 

9 

o + 



5 

NO 

10 

N + 
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The reaction system used in this report is the following: 


Reaction 

Dissociation reactions 

1 to 5 

0 2 + M= 0 + 0 + M'] 

6 to 10 

N 2 + M = N + N + M >M = 0 2 , N 2 , O, N, and : 

11 to 15 

N0 + M= N + 0 + lvJ 


Shuffle reactions 

16 

O + N 2 = NO + N 

17 

o + no=o 2 +n 


Atom-atom impact (Dissociative recombination) 

18 

0 + 0 = 0 2 + + e" 

19 

N + N = N 2 + + e" 

20 

N + O = NO+ + e" 


Electron impact (Three-body recombination) 

21 

N 2 + e~ = N2 + + e" + e" 

22 

N + e - 25 N + + e~ + e~ 

23 

O + e" — 0 + + e~ + e” 

24 

0 2 + e~ = 0 2 + + e" + e" 

25 

NO + e” = NO + + e~ + e~ 


Electron transfer (Direct neutralization) 

26 

0+0= o + + O - 

27 

O + N = N + + 0“ 

28 

O + NO = NO + + O' 


Direct ionization (Three-body recombination) 

29 

N + (O or N) = N + + e" + (O or N) 

30 

O + (O or N) = 0 + + e~ + (O or N) 
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Reaction 

Direct charge transfer for positive ions 

31 

N 2 + NO + =2 N 2 + + NO 

32 

N + NO + =2 N + + NO 

33 

O + NO + 22 . 0 + + NO 

34 

0 2 + NO + — 0 2 + + NO 

35 

n 2 + o 2 + = n 2 + + o 2 

36 

N + 0 2 + xN + + 0 2 

37 

0 + o 2 + » o + + o 2 

38 

n 2 + o + — n 2 + + O 

39 

N + 0 + — N + + O 

40 

N 2 + N + — N 2 + + N 
Electron attachment and deattachment 

41 

O’+M^O + e' + M 

42 

Og + M — 2 - Og + e + M 

43 

Og + e" 25 . 0" + O 
Rearrangement charge transfer 

44 

O + NO + 2 n + + Og 

45 

NO + NO + =2 N 2 + + Og 

46 

O + NO + — Og + + N 

47 

N + NO + =2 Ng + + O 

48 

N + NO"*" 2 0 + + N 2 

49 

NO + NO + 22 0 2 + + N 2 

50 

N 2 + 0 + 2N + + no 

51 

N0+0 + =N + + 0 2 

52 

N + 0 2 + =2 0 + + NO 

53 

O + N 2 + =2 n + + no 


Direct-charge transfer for negative ions 
54 Og + O” =2 Og” + O 
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The specific rate constants for the molecular dissociation reactions (1 to 15) have 
the following form (ref. 40): 



/SrtN 1 / 2 ia l ( D l T o 

Wl) 3» 1/2 \ T 


+ 




The specific rate constants for the remaining reactions have the form (ref. 41): 





The constants used in these equations are given in table II. 


The values for these reactions are based on information in references 40 
and 41 and were adjusted so that the theoretical rate constants would match the best 
available experimental rate data. The values were obtained from reference 42. In 

addition to references 40 and 41, some helpful information on reaction rates can be 
found in references 29 and 43 to 48. All rates are written as forward rates in the endo- 
thermic direction. Reverse rates were obtained from the following relationship: 


where the values of are computed from the free energy of the species by the use of 
a polynomial obtained by fitting the data of references 49 and 50 to a set of fourth-order 
polynomial equations for enthalpy, entropy, heat capacity, and free energy as described 
in reference 51. 


The chemical-rate terms have two forms. For reactions involving three-body 
recombination (contribution of reaction 12 to production of species 4 is a typical example), 

P^N r \ 

Q 4,12 = ^( k r,12 py 3 y 4 y 2 " k f,12 y 5 y 2) 


For the two-body reactions (contribution of reaction 20 to the production of species 6 is 
an example), 

Q 6,20 = !j"u ( k r,20 y 6 y ll ” k f,20 y 3 y 4) 

The program written to solve these equations computes and prints results for any 
desired number of streamlines in sequence. The curves in figure 12 are typical plots 
of some of the results for a body streamline at an altitude of 47.55 km. The velocity- 
altitude combinations where calculations were made are given in table m. Information 
from streamline plots such as figure 12 are used to construct profiles of gas properties. 
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APPENDIX A — Continued 
Profiles of Gas Properties in the Inviscid Flow 

Profiles of all the gas properties (p, u, T, v, and p, as well as the concentra- 
tions of the species) are found by plotting the values computed at standard network points 
against y. Initially, y values are found from the normalized plots of figure 11, but 
these values are improved by means of a computation (see appendix B) which adjusts the 
y values to ensure conservation of mass flow at the nonequilibrium conditions. 

Effect of Boundary Layer on the Electron Concentration Profiles 

Although computer programs (ref. 52) for calculating boundary- layer flow proper- 
ties along hypersonic bodies with coupled finite-rate chemistry have been developed, 
these programs are all based on the classical boundary-layer problem, for which 
boundary -layer edge properties change slowly in the flow direction and for which inviscid 
flow gradients normal to the boundary-layer edge vanish. Such idealized conditions con- 
flict with actual conditions in a thick boundary layer where vorticity due to the curved 
bow shock is large. Another deterrent to the use of coupled boundary-layer programs 
for ionized flow-field calculations arises from the mathematical complications intro- 
duced by including a large number of finite-rate reactions. 

In order to fulfill the need to compute with a reasonably complete reaction system 
and to be able at the same time to include effects of vorticity interaction, stream tubes 
from the inviscid region are continued into the boundary layer and the number of matching 
parameters along the stream tubes are increased from one to three. (This procedure is 
followed because of the necessity to delete the momentum and energy conservation equa- 
tions.) Thus, the problem reduces to determination of streamline positions in the bound- 
ary layer and of the variation of the additional matching parameters. A three-step pro- 
cedure is used to find these properties, as described in appendix C. 

Electron Concentration Profiles 

The effects of applying the procedures described herein are illustrated in figure 13. 
Curve A is a plot of inviscid stream-tube results against y-coordinates obtained from 
the mass -conservation program (appendix B). Displacement of curve A by the boundary- 
layer-displacement thickness led to curve B. The final result, curve C, includes reeval- 
uation of gas composition inside the boundary layer as well as recalculation of 
y-coordinates in the boundary layer based on merged boundary layer and inviscid pro- 
files (appendix C). 

Profile plots of N e are illustrated in figure 14. For altitudes less than 30 km, 
the only corrections made for boundary -layer effects were slight displacements in y 
of the inviscid-flow profiles. It is interesting to note that at altitudes near 40 km, the 
electron concentration is double-peaked. At higher altitudes the outer peak is dominant, 
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but at lower altitudes the peak value of N e lies very close to the body surface. Since 
collision frequency is nearly constant as a function of distance from the body, a single 
plot of v as a function of altitude for four x/djg locations is given in figure 15. Elec- 
tron concentration profiles and values of collision frequency are necessary in making 
attenuation calculations with the plane-wave program of reference 17. 

In figure 16, profile peak values of N e which include boundary-layer effects in 
the clean-air flow are plotted as a function of velocity and altitude. In the nose flow 
region, N e is essentially in local ionic equilibrium except at altitudes greater than 
70 km and at velocities less than 4 km/sec. Ionic equilibrium at low altitude is not 
maintained since velocity falls off rapidly and causes the temperature to be low which, 
in turn, means that reaction rates fall off. For local ionic equilibrium and with the addi- 
tional reasonable assumption for RAM C that y g = y + , electron concentration is 
related to atom concentration by the equation: 

N e = r e N A p = N A p ^N r O K 20 ' 

The bulge near 30 km for the aft body locations is due to local ionic equilibrium. For 
velocities less than 5.5 km/sec and altitudes greater than 40 km, the N e level on the 
aft portion of the payload will be determined by the amount of recombination of NO + + e - 
achieved in the flow time available. These two curves in figure 16 illustrate the effect 
of varying the rate constant for this reaction. The solid curve represents a reaction 
rate (ref. 48) which is three times faster than the rate given in reference 41. 

The theoretical results from references 23 and 53 are also shown in figure 16. 
Although the overall agreement is good (comparison should be made with the dashed 
lines), some mention should be made of the chemistry used in reference 23. 

The chemical reactions involving charge transfer that are used in the computations 
in reference 23 lead to errors in some of the N e values at 45 km. The only charge 
transfer reactions used were 

O + 0 2 + rz 0 + + 0 2 

N 2 + N + — N 2 + + N 

These two reactions are ineffective in removing 0 + and N + in the aft flow of the body; 
thus, the electron concentration is equal to the sum of the 0 + and N + ion concentration. 
(See ref. 1 for details.) The total positive atomic ion concentration produced at the nose, 
although not as large as the NO + concentration, does not decrease appreciably during 
expansion as does the NO + concentration. Thus, in the aft flow 0 + + N + concentrations 
are greater than that of NO + . At lower velocities NO + is the primary ion in the entire 
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flow field; thus, an incomplete charge -transfer system does not introduce any error. 
At high altitudes (>75 km) and at the velocity of 7.6 km/sec, all the reactions involving 
charged species are essentially frozen and the NO + ion is the primary ion in the entire 
flow. 
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MASS CONSERVATION COMPUTER PROGRAM 


The sketch in figure 9 should be helpful in understanding the derivation of the mass 
conservation equations. The basic idea is to equate the mass which flows through a 
specified area of the shock wave per unit time to a corresponding mass flow further down- 
stream. The mass flow inside the body of revolution generated by rotating about the body 
axis the streamline whose entry point into the shock front is distant (r s ) i from the body 
axis is pjUjir(r g ^?. 

The mass flow between the same body of revolution and the surface of the body at 
any downstream station (except those close enough to the nose to ensure that the stream- 
lines are not even approximately parallel to the body surface) is 2it * pu^r^ + y cos 0)dy. 


If these two mass flows are set equal to each other, the unknown is y^, the thick- 
ness of the gas layer. In this form the functions pu = pu(y) and ypu = ypu(y) must be 
integrated from the body surface to the streamline in question. The pitfalls inherent in 
automated curve fitting to obtain the required functions of y were avoided by formu- 
lating the problem in terms of increments of y between streamlines. For this proce- 
dure the entering mass flow is P] u i 7r |( r s)i - ( r s)i l] anc * the corresponding mass flow 


P d y * r* * 2 2 2 

downstream is 2ir J 1 pu + (y^j + y)cos 0jdy. Let A? = (r s ) i - ( r s)j_i> and 


cos 6 


+ y 


i-l* 


'0 
Then, 


Let 





. Then, 
i 


dy^_ 

2Bi 





( 2 ) 


This equation expresses the value of the increment in y, dy^, between streamlines i 
and i-l in terms of the y-coordinate for the preceding streamline y^, and the 


* 


When not subscripted, the angle 


e 


may be either 


0 b or e c . 
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integrals \]/^ and </>j. The values of p and u are already known at y = 0 and at 
y = 1, because these values are network points. By assuming that pu varies linearly 
with y in the range 0 < y < 1, the integrals can be evaluated as follows: 



The errors associated with substitution of a broken line for a smooth curve were not 
evaluated in detail but were estimated to be negligible for the purpose of calculating 
streamline position in the RAM-C flow field. 
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DETAILS OF BOUNDARY-LAYER TREATMENT 


Boundary-Layer Profiles and Displacement Effects 

Initial boundary-layer profiles are calculated from laminar similar solutions for 
an axisymmetric body at zero yaw (ref. 54) . The results from these solutions are used 
in two ways: (a) At each body station the boundary-layer displacement thickness 6* is 
added to the y-coordinate of each streamline outside the boundary layer (6* is approxi- 
mate since the computation does not include vorticity) , and (b) the profiles of pu and u 
are merged with corresponding profiles in the displaced inviscid-flow field. Nonequilib- 
rium edge conditions are used, and the program assumes that boundary-layer gas proper- 
ties are thermally frozen (that is, no dissociation or recombination of atoms occurs in 
the boundary layer) . This assumption is reasonably good for finding thermal properties 
(not ionization), as chemical kinetics computations show that very little change in atom 
concentration occurs in most of the boundary layer aft of the shoulder in the altitude 
regime where boundary- layer effects are significant. Even near the wall where atom 
recombination will actually occur, the error is very small, since T and p in this 
region must approach the cold-wall values as specified by the ablation temperature. Use 
of local nonequilibrium edge properties for the boundary layer gives better results than 
use of either equilibrium or frozen-flow edge properties (as in a complete blunt-body 
program), since, in general, processes behind the shock and in the expansion around the 
body are far from either. 


Initial nonequilibrium edge properties are taken from the inviscid body streamline 
calculations. Then, since the boundary layer is thick, swallowing of the inviscid stream- 
lines into the boundary layer is taken into account in choosing edge conditions along the 
boundary layer for a repeat calculation. The criterion for swallowing a given streamline 
is based on equality of the boundary-layer mass flow to the entering mass flow. For 
stations in the nose region, let y T = y^l + Then 


Pl u l ir ( r s)i = ffd N cos e b J 0 i° u <*y’ 

( cos On v \ 

1+ . Then 

2 r b / 


( 5 ) 



( 6 ) 
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The first use of the boundary- layer solutions is to provide values of 6* which are added 
to the body dimensions and thus displace the inviscid shock layer. A plot of the new shock 
shape is used to provide new starting angles and pressure for the entering streamlines. 
Also, an effective cone angle and nose diameter are determined which, along with use of 
the normalized flow-field inputs previously discussed, allows new pressure distributions 
along the streamlines to be found. The new pressure distributions are used to recompute 
the inviscid nonequilibrium shock-layer properties and streamline locations, and the 
whole cycle can be repeated until satisfactory agreement is obtained between input and 
output shock shape and streamline position. 

Merged Profiles 

Flow vorticity was accounted for partially in the preceding step, where edge proper- 
ties were chosen by means of a swallowing criterion. An additional correction for the 
vorticity normal to the flow direction is obtained by merging boundary-layer profiles with 
displaced inviscid-flow profiles. 

Merging is carried out for the pu and u shock-layer profiles, and three princi- 
pal criteria are applied in the merging procedure: (a) mass conservation for the pu 
profiles, (b) velocity defect for the u profiles, and (c) smooth transition for both from 
the displaced inviscid profile (in the neighborhood of y = 6) to the boundary -layer profile 
(in the neighborhood of y = 6*). In figure 17(a) where a typical pu profile has been 
plotted against the transformed coordinate y* , equality of the shaded areas indicates that 
the mass flow is the same whether obtained by integrating the merged or the displaced 
inviscid profile. In addition to respecting the mass-flow criterion, the merged profile 
is required to blend smoothly into the two original profiles. 

Figure 17(b) illustrates the merging of the velocity profiles. The velocity defect 
inside 5 m (s m was determined in step (a)) on the merged profile is made similar to 
that within the boundary-layer profile. This velocity defect criterion is applied only for 
the outer region of the profile, since the inner profile must merge into the boundary-layer 
profile. 


Boundary-Layer Stream Tubes and Matching Conditions 

The next step in the flow-field determination is to identify "streamline" positions 
inside the merged boundary layer so that the matching parameters for use in the boundary- 
layer nonequilibrium stream -tube analysis may be evaluated. This evaluation is made by 
integrating the pu against y’ merged profiles to find those values of y' within which 
mass flow equals that of each entering streamline. With streamline positions established 
in the boundary layer, values of the matching conditions are read from merged plots of 
pu and u at each body station. These values, along with the pressure-matching 
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parameters (which are assumed to be unchanged from those of the corresponding dis- 
placed inviscid streamlines), are used in the boundary- layer version of the stream -tube 
program. Profiles of species concentration can then be constructed at each body station. 
It should be pointed out that the specification of p and p in the boundary- layer stream- 
tube analysis is nearly equivalent to specifying the temperature T, within the altitude 
range where the boundary layer is significant. In this range the composition changes in 
the boundary layer involve principally electrons and ions, with no significant alteration 
of the neutral species concentration. Because the chemical kinetics depends strongly on 
temperature, it is important that the boundary layer and matching conditions be deter- 
mined by using finite-rate considerations so that the temperature in the boundary layer 
is correctly specified. 

The profiles of electron and ion concentrations which are obtained by these methods 
do not contain effects of surface recombination or of diffusion of ionic species, and there- 
fore restrict the validity of calculations for the body considered herein to altitudes below 
about 70 km. 



APPENDIX D 


EFFECTS OF ABLATION ON ELECTRON CONCENTRATION 
IN BOUNDARY LAYER 

Theoretical estimates are made of the influence on electron concentration of the 
ionization of NARMCO 4028 ablation products in the RAM C-I boundary layer. These 
estimates are based on a calculated rate of ablation, on measurements of the ionizable 
impurities present in NARMCO, on theoretical distributions of air flow and ablation 
material in the boundary layer, and on assumed local equilibrium ionization of the abla- 
tion impurities. 


Ablation Gas Flow Rate 

The rate at which the surface of the NARMCO nose cap material recedes due to 
ablation, and the flux of the vaporized and pyrolyzed material into the boundary layer at 
various stations on the body, has been computed by Ernest V. Zoby of Langley Research 
Center over the complete entry altitude range for RAM-C. These data were then cross- 
plotted, at various given altitudes, in the form of m, the ablating surface flux, against A, 
the integrated surface area of the nose cap to the several body stations, 

A = A(x) 

The faired plots were then graphically integrated 


m 


abl 


= f 


m dA 


(V 


to find the ablation flow rate m abl past the body stations x/d^j = 0.15, 0.422, and 0.81 
at each of the given altitudes. The resulting values of m^ are plotted as a function 
of altitude in figure 18. The station x/dj^ = 0.81 represents the termination of the 
NARMCO nose cap and the beginning of the teflon body coating. Since only the NARMCO 
is herein considered as a source of ionizable material, the values for m a j^ at 
X /^N = 0*^1 are also used for all other downstream body locations. It is seen in fig- 
ure 18 that the peak ablation rate at x/d^ = 0.81 occurs at an altitude of about 25 km. 


Ablation Gas Flow Distribution 

To determine the distribution of ablation pyrolysis gas flow in the RAM-C flow 
field, it is first assumed that the ablation gases are confined to the boundary layer, and 
that the amount of the gases approaches zero at the outer edge. This assumption means 
that integration of the ablation flow across the boundary layer at a given station must be 
equal to the value of m a ^ already determined for that station (see sketch (a)). 
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abl 


or 


r* y=5 

m abi = J =0 (p u ) ab i M R b + y cos e c) d y 


m abl = 27rR bI 0 6 f 1 + 


^cos 0 c )(pu) abl dy 


This expression can be written in terms of a transformed normal coordinate y' as 

r 


m abl = 27rR bJ 0 (P u )abl d y’ 


( 8 ) 


( 9 ) 


where 


y ' = y ( 1 + 2Rg COS,, c) 


dy' = dy f 1 + -¥~ cos 6 C 


where the ordinate is proportional to the flow cross-sectional area. If (pu) ab ^ (see 
sketch (b)) is defined as 
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y’ 



~ (P u )abl 



and plots of (pu/pu) abl against y’/6’ , the normalized distribution of ablation gas flow, 
must have an integrated area of 1.0. 

The second principal assumption is that the form of the normalized distribution be 
generally similar to the form of the velocity defect in the air boundary layer (l - jj-) 

/ x \ ^ 

plotted against y'/6' for stations along the nose cap ^ 0.810J and that for stations 

downstream of the nose cap, there is a peak in the distribution which moves progressively 
out from the surface toward the point of maximum shear, the pu value falling off rapidly 
at the surface as in reference 55. The latter characteristic reflects the fact that only the 
NARMCO is considered to be pertinent to the ablation ionization problem. The assumed 
form which is used herein is shown in figure 19. 

By using values of the boundary-layer thickness, 

5’ = 6(1 + — cos 0 C ) 

\ / 

and m abl , both at given body stations, the value of (pu) abl is found from equation (10) 
and values of (pu) abl plotted against y can then be found by using figure 18. 
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Air Boundary-Layer Profiles 

The distributions in the air boundary layer of the parameters pu and u are 
found by use of a merging technique described in appendix C which partially accounts for 
vorticity interaction of the boundary layer with the inviscid nonuniform flow field. The 
technique merges the displaced inviscid-flow profiles with the classical laminar boundary- 
layer profiles in such a way as to provide a smooth transition from the viscous flow in 
the inner part of the boundary layer to the inviscid flow at the outer edge of the boundary 
layer. 

The temperature profile in the low-altitude boundary layer T = T(y) is computed 
by using the merged pu and u values at a given y-value 

T £i! (12) 

Z i^ U) 

where p, the pressure, is from the RAM-C shock layer solutions, and Z, the compres- 
sibility factor, is assumed to vary smoothly (from faired plots) from the nonequilibrium 
values at the outer part of the boundary layer to a value of 1.0 at the wall (y = 0). This 
method to obtain temperature was used at low altitudes where there are not enough stream- 
lines in the boundary layer to define an accurate temperature profile. It is seen from 
equation (12) that the value of T becomes zero or indeterminant near the wall for body 
stations greater than x/d^f = 0.81. This result presents no particular problem, however, 
because the wall is assumed to be at 1000 K and T is simply faired toward this value. 

The arbitrariness is not important since no significant ionization of the ablation impurities 
occurs at the lower temperature (T S 1500 K). 

The use of laminar boundary-layer theory in the merging technique is believed to 
be generally appropriate for these ablation calculations. The transition from laminar to 
turbulent flow at the x/d^ = 0.81 station is not expected to start above an altitude of 
about 28 km. At lower altitudes, only the thickness and distribution of ablation is in 
error, since the ablation rate calculations included the effects of turbulent flow. The 
resulting error in calculated radio signal attenuation is not known but may be small 
because of the compensating effects of larger thickness but lower (pu) ^ values (lower 
N e ) for turbulent flow. (See eq. (10).) 

Ionization of Ablation Gas in Boundary Layer 

In order to determine meaningful estimates of the ionization within reasonable 
limitations on computing effort, several simplifying assumptions are employed. It is 
first assumed that the ionization is due solely to the alkali impurities known to be present 
in NARMCO. This assumption is probably good, since these impurities are present in 
amounts at least to 10“3 ? by weight, and the ionization potential for the alkalis is very 
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much lower than that of other species expected to be present. Reference 56 gives the 
species that would be expected in the boundary layer when a phenolic carbon heat shield 
ablates. The second assumption is that the pyrolysis gases injected into the boundary 
layer contain these alkali impurities to the same extent as the virgin NARMCO material 
which was tested. The validity of this assumption is not known. The third assumption 
is that the ionization of the alkali impurities is in local equilibrium with the nonequilibrium 
air-ablation gas mixture in the boundary layer. This last assumption would be expected 
to be reasonably good for blunt-nosed bodies at low altitudes where large numbers of gas 
collisions occur in the flow field. However, at intermediate and high altitudes, the local 
equilibrium assumption would become progressively worse since there would not be time 
for recombination of electrons and ions in the expanded flow at intermediate altitudes, 
and, at the higher altitudes, there would be a lack of equilibrium ionization in the nose 
region as well. On the other hand, theoretical nonequilibrium calculations (ref. 56) for 
the ablation-air mixture would be subject to great uncertainty due to lack of knowledge 
with respect to the reaction rates for the many possible reactions. Because of the large 
computing effort required for nonequilibrium calculations, and the other uncertainties 
involved in the boundary-layer calculations, the local-equilibrium assumption is used for 
simplicity, its limitations being recognized for later quantitative comparison with flight 
data. 


The three flow parameters which are needed for the calculations are 


<P u )abl 

(pu) 


air 


p/p Q , and T. These parameters have already been discussed, but it should be pointed 

(P u ) a bl 

becomes indeterminate at y — 0, for the body stations beyond 


out that 


(pu) 


air 


x /% = 0.81. This problem is only academic because near the wall the ionization is well 
below the peak and is falling rapidly toward very low values. 


The Saha equation is used to compute the ionization of the boundary- layer mixture, 


qj2 

1 - 0 ! 



log 


-1 


2.5 log T - 5040 


6.49 + log G 


l 


X- — 

A 1 p 


(13) 


where G, the degeneracy factor, is assumed to be 1 and hence the term drops out. A 
value of 4.85 eV is used for Vj, the ionization potential, and represents a value between 


that for pure Na and pure K. 


The parameter 



ionizable (alkali) species and is computed from 


represents the mole fraction of the 
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where 


and 


I 


X A = X 


V X a ik 

abl ^ x aW 


3*. 

Xair ^ u >air\ M abl/ 

1 + ^abl ^ + (p u ) abl /Mair\ 

x air (pu) a ir yMabl/ 



M abl\ 

M alk/i 


(14) 


(15) 


(16) 


Combining equations (14), (15), and (16) and substituting the values for the various molec- 
ular weights gives 


0 ? ( ^ u) abl V 

Y x 

^ 1 " (P^)abl 


(17) 


1 + 3.5 


(pu) 


air 


where M al ^ = 30, a value between that of Na and K, was used; M a i r = 21, a value corre- 
sponding to an overall average for the air boundary-layer flow; M a bi = 6, a value arbi- 
trarily assumed since the composition of the NARMCO was not known. 

V (P u ) a bl 

By using values of > Fi and values of - — r plotted against y, the values of 

*-• ( pu) air 

^ X^ plotted against y can be obtained by use of equation (17). Then a plotted 

against y can be obtained by use of equation (13). The electron mole fraction and elec- 
tron concentration are then obtained, respectively, from 


X, 


-I*. 


N e = 


7.35 X 10 21 ( 

w 


The resulting computations of N e plotted against y for the ablation impurity 
ionization in the RAM-C boundary layer show large effects due to the NARMCO ablation, 
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as is seen by the examples shown in figure 20. Diffusion effects were not included in 
obtaining the N e profiles since they are not believed to influence N g p significantly 
for these nondissipative blunt -body boundary layers at low altitude. Calculations were 
made for x/d^ locations which bracket RAM C-I flight data. 
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INTERPRETATION OF LOW- ALTITUDE REFLECTOMETER DATA 

It was pointed out in reference 8 that determination of the flight time at which the 
peak N e at a given body station goes through the critical value is made by using the 
sudden rise or decay of reflection coefficient. Although this concept was shown to be 
straightforward for the high-altitude data, it was pointed out that for low altitudes there 
were uncertainties due to high collision frequency, rapid erosion of the heat shield, body 
motion, antenna detuning, and so forth. Several of these factors are discussed further. 

Heat-Shield Erosion 

In figure 3, which shows the reflection coefficient data, as given in reference 8, the 
regions of fast ablation are indicated. Some ablation occurs prior to and subsequent to 
these times but the rates are much lower. The total changes in teflon thickness occurring 
during the times shown for body stations 1 to 4 are 1.9, 0.47, 0.29, and 0.23 cm, respec- 
tively, according to calculations made by Ernest V. Zoby of Langley Research Center. 

The free-space wavelengths for the K-band, X-band, S-band, and L-band reflectometers 
are 0.86, 3.0, 9.0, and 27 cm, respectively. Effects due to thickness change for the Kq 
reflectometer are shown in figure 21 where the calculated rate of change of thickness, 
normalized to the peak rate, is compared with the normalized observed frequency of the 
Kq Tp oscillations which are seen in figure 3(a) . 

The close correspondence between the two parameters shows that the reflection 
coefficient for this reflectometer depended strongly on teflon thickness, since several 
swings of Tp from values near 1.0 to values near 0 are seen to occur during the abla- 
tion period but within a period when the plasma is maintained far above the critical N e . 
Note that the Tp swings occurring prior to the fast ablation period are for a calculated 
change of only 0.18 cm in ablation thickness, during which time the plasma is well above 
critical N e . It is also interesting to note that the magnitude of the Ip swings decreases 
to a relatively small value at about the time when the plasma N e p falls through the 
critical value according to the plasma theory. There are similar swings, or bumps, but 
of lower magnitude and frequency, seen on some of the other records during the fast 
ablation period. For example, the bumps in Xj_, Sj, and Xg superimposed on the decay 
region of the records make determination of the actual decay very uncertain. One could 
sketch in any number of possible shapes which might exist in the absence of ablation, as 
shown by the dashed shapes inserted in figures 3(a) and 3(b) . 
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Body Motions 

There was a brief period during the low-altitude part of the RAM C-II reentry that 
the vehicle went from an angle of attack of around 1° to values as high as 5°. This body- 
motion period is seen in figure 22 and is from about 414.0 to 418.5 seconds, or from an 
altitude of about 38 to 30 km. In order to make estimates of the changes in N e during 
this period, an approximate flow-field theory was employed. For reflectometer sta- 
tions 2,3, and 4, it was assumed that the flow conditions along the most windward and 
most leeward rays of the conical body at angle of attack could be approximated by com- 
puting axisymmetric flow fields for 14° and 4° half-angle sphere cones, respectively, to 
examine the case of an angle of attack of 5°. For station 1, conditions at angle of attack 
were approximated from the zero angle case by computing for stations 5° forward and 
aft of the actual body station angle (angle between body axis and line from sphere center 
to body station) . The computed peak electron concentrations at each body station for 
altitudes of 30.7 and 40 km, which approximately bracket the body-motion region, are 
shown in figure 23. The hatched area is the resulting estimated swing in N 6j p between 
the windward (upper) and leeward (lower) limits due to spin at an angle of attack of 5° 
during the body-motion period. The N e profiles for the 14° and 4° cases are shown in 
figure 24. 

It can be seen that for some of the body stations, a large change in peak electron 
concentration is predicted at angle of attack as the measurement point on the body goes 
from leeward to windward during each body rotation (spin rate of several times per 
second). Also the profile shape changes from leeward to windward. For the station 1 
case in figure 23, the N e p change is not so large and this, along with the fact that the 
N e level is far above the critical values for the reflectometers, would suggest little, if 
any, oscillation of the reflection coefficient records due to N e change during the body- 
motion period. Little change was observed. 

For station 2, a fairly large oscillation of Tp is suggested for K-band, since N e 
swings through critical, and this oscillation is indeed seen (fig. 23). However, changes 
in N e profile shape can also produce changes in T p and hence cannot be separated 
from the N e p effects. The other frequencies show a small but discernible effect 
during the body-motion period as might be expected since the N e values remain sub- 
stantially higher than critical. However, for stations 3 and 4 there is a very large swing 
of theoretical peak N e due to body motion and for X-band, in particular, where the 
swing is through critical N e , very large effects on reflection coefficient were observed. 
Note that the hatched region also goes above K-band critical for about 2 seconds near the 
end of the period, but it must be recalled that the hatched band corresponds to an angle 
of attack of 5° throughout the body-motion period. As seen from figure 22, the angle is 
below 5° during this 2-second time so that K-band effects due to N e?p are not really 


38 



APPENDIX E - Continued 


predicted. S-band effects are somewhat smaller since the swing is not through critical, 
and since the spatial changes in profile shape are a smaller fraction of a wavelength. 

Boundary-Layer Transition 

Careful inspection of the reflection-coefficient records in figure 3 shows that in 
each of the four X-band records, one very sharp decay of reflection coefficient occurs 
which is qualitatively unlike the decays in the other instruments. (The K 2 record is 
discussed later.) That is, these decays are very sharp, nearly discontinuous, even on 
the time -amplified records given in reference 8 . The times at which these sudden decays 
occur for stations 1 to 4 are 420.5, 419.3, 418.1, and 415.6 seconds, respectively. The 
sudden decay on X 2 is only different from the other three in that the Tp does not drop 
all the way to a near-zero value. Another way of noting the uniqueness, or suddenness, 
of these decays is to compare them with the slopes of the S-band and L-band decays, or — 
more dramatically — with the slope of the X 2 final decay, which is significantly less steep. 
(See time -amplified records of ref. 8 .) 

It is interesting to compare the altitudes at which these sudden decays occur with 
those altitudes at which boundary-layer transition would occur based on a local Reynolds 
number criterion of 350 000. This value is the value assumed by Zoby in making esti- 
mates (believed to be conservative) of the RAM-C heating rates. The comparison is 
shown in figure 25 where it can be seen that very good correlation is obtained. This cor- 
relation is believed to be a strong indication that these sudden decays of Ip are asso- 
ciated with transition. They are not associated with the body motions or fast ablation 
periods since they occur at progressively different times with respect to these periods, 
and none occurs exactly at the beginning or end of the periods. The first two decays (X 4 
and X 3 ) occur during the body-motion period and may be somewhat different in charac- 
ter than the last two decays since the instrument rotates a few times in and out of transi- 
tion when spinning at an angle of attack (because of different flow lengths from stagnation 
point). 

Since the important influence which heat-shield erosion, body motions, and boundary- 
layer transition can have on the low-altitude reflectometer records has been seen, there 
remain some questions to be clarified. One question is in regard to the sudden decay of 
the K 2 reflection coefficient, which occurs neither at the transition time nor very near to 
the time predicted for critical N e from the theory. (Note the arrows in fig. 3.) How- 
ever, it does occur at a time coincident with the end of the body-motion period. It is 
believed that the decay might be qualitatively explained in terms of body motion and heat- 
shield erosion as follows: (a) It was previously shown that the Kj reflection coefficient 
went down markedly when the surface ablation was only on the order of 0.1 cm (that is, 
first drop occurring before fast ablation), (b) During the fast-ablation period at 
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station 2, there is theoretically a thickness change of about 0.47 cm and it is reasonable 
that approximately 0.1 cm occurs from the start of fast ablation to the end of body motions. 
(See fig. 3(b).) (c) If, therefore, a large r p oscillation due to body motion is super- 

imposed on a decaying r p due to ablation, it is possible that the r p indication could 
cease at a time coincident with the end of large body motion. It is seen from K 2 in fig- 
ures 3(b) and 23(a) that the first r p rise (which occurs at 414.7 sec) is almost exactly 
at the time predicted from theory which includes body motion. Although the magnitude 
of the oscillations should probably then increase, because of the general increase in 
N e p , any drop in r p due to ablation would be working against this increase. At the 
time when body motion ceases, the r p reduction from ablation may already be quite 
large. However, this argument is only qualitative based on the assumption that K 2 is 
generally similar to Kj in its sensitivity to heat-shield thickness. 

There is an important question with regard to the sudden drops in T p which are 
believed to be associated with transition, since both changes in N e p and/or profile 
shape can produce changes in reflection coefficient. A drop in N e>p through critical 
would certainly result in a drop of T p to a very low value, whereas a shift in the N e 
profile shape could cause r p to go through a null point but it would probably not drop 
to a low value for those cases where the reflectometer is extended to the surface. (See 
fig. 8 in ref. 8.) Yet there is no conceivable way that the N e p can fall to near critical, 
for example, at transition for station 1. Before transition the laminar boundary-layer 
thickness at this location is theoretically about 0.15 cm, and the plasma is nearly three 
decades supercritical for the X-band all the way out to the shock wave (about 1.7 cm) as 
seen in figure 14. Even if the boundary layer was several times thicker after transition 
and if it is assumed that teflon in the boundary layer is an efficient electrophilic, the peak 
value of N e in the inviscid plasma is virtually unchanged. Furthermore, a thickening 
of the boundary layer due to blowing (the blowing parameter is theoretically about 0.06) 
or surface roughness would simply move the shock and inviscid plasma further out, with 
negligible effect on N e>p . 

There is an additional point which suggests that the large sudden drop in the Xi 
record is more likely due to a profile change than to a large drop in N ep . A large drop 
in N e p in the region around X-band critical should presumably also show up as a sig- 
nificant drop in the record since the N e change was at a level of about 10 times 
S-band critical and collision frequency was high (that is, r p was not at its highest 
value). For example, during the body-motion period, there were changes in the S-band 
records when N e was supercritical for S-band and swinging about the X-band critical 
level. These S-band changes were probably not due to profile shifts, because the 
boundary-layer changes at low altitude are small compared with the S-band wavelength. 

By the same reasoning, little change would be expected in Sj if a sudden change occurred 
in the distance of the N e profile from the body surface but without any change in N e jP - 
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(This change corresponds to a sudden thickening of the boundary layer for transition from 
laminar to turbulent flow.) Thus it seems significant that S-band response is missing at 
the times designated for transition, even though X-band response is large, because the 
changes in boundary-layer thickness are not negligible compared with the X-band 
wavelength. 


High Electron Collision Frequency 

In figure 26 the theoretical effect of collision frequency on the variation of Tp 
with N e /N e cr is shown for the X-band reflectometer used at station 1 on RAM C-n. 
(Similar curves for the X-band reflectometers used at stations 2,3, and 4 are given in 
ref. 8.) These calculations are based on a parabolic N e profile where the 4n p n 
point is located about 1/4 cm from the surface of the vehicle. Although the actual flight 
profiles are not known, the theoretical profiles shown in figure 14 (including boundary 
layer) for the low altitudes pertinent to this discussion have N e> p points which are in 
this general range. In figure 15 it can be seen that the theoretical value of v for sta- 
tion 1 at the altitude for theoretical N 6jCr (approximately 20 km) is about 4 x lO-^- sec - ^. 
According to figure 26, the value of Tp at N e = N e cr is less than 0.01 and at 
N e = 10N e>cr is only about 0.1. It is obvious that identification of N e>cr from a record 
of low-altitude decay of Tp for this case would be very uncertain, even within wide lim- 
its for N e> p. More exact calculations, based on the theoretical plasma profiles, could be 
made but the large effort would not be warranted merely to show qualitatively the strong 
effect of collision frequency in the records as discussed here. 

It would be interesting if a calculation of Tp based on the theoretical plasma pro- 
files and reflectometer characteristics could be made for comparison with all the records. 
Although it would be possible to make these calculations for the X-band reflectometers, 
with a considerable effort, they have not been made. As pointed out previously, the pro- 
files used to compute the X-band characteristics in figure 26 and reference 8 are not 
grossly different from the theoretical plasma profiles at low altitude where the decays 
occur. Therefore, for the purpose of making at least qualitative comparisons in the 
decay region for the X-band, values of Tp were calculated by using reference 8 and fig- 
ures 2, 26, and 15. The results are shown in figure 27(a), where the numbers beside the 
squares are values of Ngy'Ng cr . It is seen that the only station where the comparison 
is at all close is at station 2 where the transition apparently did not cause a complete 
decay of Tp as at the other stations. 

In order to make similar, but even cruder, comparisons with the S-band and L-band 
records, a very approximate approach was used. For this purpose the high collision fre- 
quency curves in figure 26 and reference 8 for X-band were used in normalized form 
(r p plotted against v/w and N e /N e>CI A for application at the lower frequencies. The 
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results are shown in figures 27(b) and 27(c). In spite of the admittedly crude approach 
used, it is believed that the results serve to illustrate the order of uncertainty in inter- 
pretation of the low-altitude results. That is, at low altitudes where v/w > 1, the steep 
drops in reflection coefficient appear to be associated with changes in N e /N e cr between 
100 and 10, and between 10 and 1. Thus, the region of steep drop in r p at low altitude 
does not necessarily include the N e = N e>cr point, because of the effects of high colli- 
sion frequency. 

In view of the discussions in this appendix, the reflectometer decay bars in figure 2 
are shown dashed to indicate these uncertainties, in contrast to the solid rise bars at 
higher altitudes which are believed to be much more precise indicators of the N e = N e cr 
point. 
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TABLE I.- LOCATION OF ONBOARD RADIO FREQUENCY AND ELECTROSTATIC PROBES 


x/d N 

(a) 


0.15 (45°) 


RAM C-I 


(Station 1) 
0.76 


(Station 2) 
2J3 


(Station 3) 
2.4 

2.6 


VHF telemetry slots 
(259.7 MHz) 
X-band telemetry 
horns (9210 MHz) 


3.2 


RAM C-II 

K-band reflectometer (35 000 MHz) 
X-band reflectometer (10 044 MHz) 
S-band reflectometer (3 348 MHz) 

K-band reflectometer (35 000 MHz) 
X-band reflectometer (10 044 MHz) 
S-band reflectometer (3 348 MHz) 
L-band reflectometer (1 116 MHz) 

K-band reflectometer (35 000 MHz) 
X-band reflectometer (10 044 MHz) 
S-band reflectometer (3 348 MHz) 
L-band reflectometer (1 116 MHz) 


X-band telemetry horns (9 210 MHz) 
VHF telemetry ring (259.7 MHz) 


3.4 


3.5 


VHF telemetry ring 
(225.7 MHz) 


(Station 4) 

3.5 

3.7 C-band beacon horn 

(5700 MHz) 


K-band reflectometer (35 000 MHz) 
X-band reflectometer (10 044 MHz) 
S-band reflectometer (3 348 MHz) 
L-band reflectometer (1 116 MHz) 

C-band beacon horn (5 800 MHz) 


3.8 

3.9 


Langmuir probe 
rake* 3 


VHF telemetry ring (225.7 MHz) 


4.0 


Langmuir probe rake 


b 


a Location of center line of component, 

d N = 31.90 cm (12.56 in.) for RAM C-I NARMCO 
= 30.48 cm (12.00 in.) for RAM C-H beryllium 
= 32.13 cm (12.65 in.) for RAM C-H teflon. 
b Probes located at y = 0.96, 1.84, 2.70, 3.58, 4.46, 5.34, 6.20, and 7.08 cm. 
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TABLE II.- CONSTANTS 


l 

Reaction 

H 


p 

0 ^ , cm^ 

1 

C>2 + C>2 — 0 + O + 0 2 

16.000 

217.34 

14.60 X 10- 

2 

02 + N 2 = 0+ 0 + N 2 

14.934 

217.34 

7.75 

3 

O 2 + 0 — ^ 0 + 0 + 0 

10.667 

217.34 

60.00 

4 

0 2 +N s o+o + n 

9.739 

217.34 

6.72 

5 

0 2 + NO 0 + 0 + NO 

15.484 

217.34 

5.00 

6 

N 2 + 0 2 2 = N + N + 0 2 

14.934 

414.50 

10.30 

7 

N 2 + N 2 ss N + N + N 2 

14.000 

414.50 

16.50 

8 

N 2 +0sN+N+0 

10.182 

414.50 

6.45 

9 

n 2 +n~n + n + n 

9.334 

414.50 

72.00 

10 

N 2 + NO N + N + NO 

14.483 

414.50 

7.00 

11 

NO + 0 2 — N + 0 + 0 2 

15.484 

275.74 

2.94 

12 

NO + N 2 N + 0 + N 2 

14.483 

275.74 

3.12 

13 

NO + 0 ^ N + 0+ 0 

10.435 

275.74 

6.06 

14 

NO + N ss. n+O + N 

9.545 

275.74 

11.30 

15 

NO + NO — N + 0 + NO 

15.000 

275.74 

14.80 

16 

0 + N 2 -N0+N 

10.182 

138.93 

2.69 

17 

O + NO =■ 0 2 + N 

10.435 

71.24 

1.10 

18 

0 + 0 0 2 + e 

8.000 

294.85 

.175 

19 

N + N N 2 + + e“ 

7.000 

247.27 

.66 

20 

N + 0 2= N0 + + e" 

7.467 

117.69 

.076 

21 

N 2 + e =• Ng + + e" + e™ 

5.4858 X 1G“ 4 

661.90 

1.00 

22 

N + e“ ~ N + + e' + e“ 

5.4858 X 10- 4 

617.80 

1.00 

23 

0 + e’ = 0 + + e~ + e" 

5.4858 X 10- 4 

578.40 

1.00 

24 

0 2 + e ^ 0 2 + e + e 

5.4858 X 10~ 4 

513.00 

1.00 

25 

NO + e~ ^ N0 + + e” + e" 

5.4858 X 10" 4 

393.00 

1.00 

26 

0 + 0 =■ o + + O' 

8.000 

516.00 

1.00 

27 

0 + N — N + + O' 

7.467 

555.50 

1.00 


FORMULAS 


l 

Reaction 



o-j , cm 25 

28 

O + NO — NO + + O" 

10.435 

330.70 

1.00 x 10“ 16 

29 

N + (O or N) — N + + e" + (O or N) 

7.000 

617.80 

1.00 

30 

O + (O or N) — 0 + + e" + (O or N) 

8.000 

578.40 

1.00 

31 

N 2 + NO + =s N 2 + + NO 

14.50 

268.90 

20.00 

32 

N + NO + — N + + NO 

9.55 

224.80 

20.00 

33 

O + NO + = 0 + + NO 

10.40 

185.2 

20.00 

34 

0 2 + NO + = 0 2 + + NO 

15.50 

120.20 

20.00 

35 

n 2+ o 2 + = n 2 + + o 2 

14.90 

148.70 

20.00 

36 

N + 0 2 + = N + + 0 2 

9.74 

104.50 

20.00 

37 

0 + o 2 + = o + + o 2 

10.70 

65.00 

20.00 

38 

N 2 + o + = n 2 + + 0 

10.20 

83.70 

20.00 

39 

N + 0 + = N + + O 

7.47 

39.51 

20.00 

40 

N 2 + N + = N 2 + + N 

9.34 

44.19 

20.00 

41 

0” + M — 0 + e"+M 

10.00 

62.48 

10.00 

42 

0 2 ‘ + M 0 2 + e“ + M 

15.00 

19.55 

10.00 

43 

0 2 + e~ 2 =. O' + O 

5.4858 X 10" 4 

155.14 

.10 

44 

O + NO + ss N + + O 2 

10.44 

282.90 

20.00 

45 

NO + NO + =■ N 2 + + 0 2 

15.00 

188.60 

20.00 

46 

O + NO + — 0 2 + + N 

10.40 

178.40 

.20.00 

47 

N + NO + — N 2 + + O 

9.55 

130.00 

20.00 

48 

N + NO + =. o + + N 2 

9.55 

46.31 

20.00 

49 

NO + NO + = 0 2 + + N 2 

15.00 

39.93 

20.00 

50 

N 2 + 0 + = N + + NO 

10.18 

178.40 

20.00 

51 

NO + 0 + — N + + 0 2 

10.44 

97.71 

20.00 

52 

N + 0 2 + =. 0 + + NO 

9.74 

6.80 

20.00 

53 

O + N 2 + = N + + NO 

10.18 

94.74 

20.00 

54 I 

O n + O — ► Oo + O 

10.67 

42.93 

20.00 



TABLE III.- 


Time from 
lift-off, sec 


426.8 

19.202 

424.8 

21.336 

423.3 

22.860 

422.3 

23.774 

421.4 

25.146 

419.4 

27.999 

418.0 

30.632 

416.4 

33.706 

413.1 

40.000 

409.2 

47.549 

406.2 

53.340 

402.0 

61.570 

397.3 

71.018 

392.6 

80.162 

390.2 

85.040 


C TRAJECTORY DATA 


ude 

Velocity 

ft 

km/sec 

ft/sec 

63 000 

3.048 

10 000 

70 000 

3.883 

12 740 

75 000 

4.572 

15 000 

78 000 

5.029 

16 500 

82 500 

5.425 

17 800 

91 860 

6.218 

20 400 

100 500 

6.614 

21 700 

110 583 

6.980 

22 900 

131 233 

7.391 

24 250 

156 000 

7.565 

24 820 

175 000 

7.620 

25 000 

202 000 

7.651 

25 100 

233 000 

7.651 

25 100 

263 000 

7.651 

25 100 

279 000 

7.651 

25 100 



Time, sec 



i I I L 

4 5 6 7 

Velocity, km/sec 

Figure 1.- RAM C configurations and entry profile. 
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Figure 2.- Comparison of peak electron concentration (air theory) 
with RAM C-II flight results. 
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Power reflection coefficient 



(a) Station 1 = 0.15). 

Figure 3.- Measured reflection coefficients taken from reference 8. 



Power reflection coefficient, 


cr 



(b) Station 2 = 0.76^. 

Figure 3.- Continued. 
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Power reflection coefficient. 
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(c) Station 3 ss 2.3^. 

Figure 3.- Continued. 
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Power reflection coefficient, 


K -band 



(d) Station 4 = 3.5^. 

Figure 3.- Concluded. 
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Phase of reflected wave, deg 
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(a) Station 1 (x/djj = 0.15). 

Figure 4.- Measured phase coefficients taken from reference 8. 
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Phase of reflected wave, deg 


K -band 



Time from lift-off, sec 


(b) Station 2 (x/djq- = 0.76^. 
Figure 4.- Continued. 
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Phase of reflected wave, deg 


K -band 



(c) Station 3 (x/d^r = 2. 3 ). 
Figure 4.- Continued. 
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Phase of reflected wave, deg 
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(d) Station 4 (x/djq- = 3.5). 
Figure 4.- Concluded. 
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Signal loss, dB 



Figure 8.- VHF attenuation (226 MHz) for x/d-^f = 4 using plane wave 
theory for RAM C-II vehicle at angles of attack of 0° and 5°. 
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Figure 13.- Steps in calculation of an electron concen- 
tration profile. Altitude, 71.02 km; x/d^ = 4. 



y. cm 


(a) Altitude, 21.336 km. (b) Altitude, 23.774 km. 

Figure 14.- Electron concentration profiles for air. (Fast NO + + e 
recombination rate of ref. 48 used.) 






(g) Altitude, 40.000 km. 


(h) Altitude, 47.549 km. 
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Figure 14.- Continued. 
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Figure 16.- Comparison of computed peak electron concentration 
(air theory) with previous calculations. 
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Figure 22.- Payload angle-of-attack time history. 
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Figure 25.- Correlation of reflectometer decay 
from RAM C-II flight with boundary-layer 
transition. 
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Figure 26.- Theoretical variation of reflection coefficient for 
X-band reflectometer at station 1. 
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(a) X-band. 


Figure 27.- Comparison of estimated reflection coefficient 
(Numbers beside estimated points are values of 


decay with experiment. 
^e/^e,cr-) 
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Experiment 

---D--- Faired estimated values 






(b) S-band. 

Figure 27.- Continued. 
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-o- Faired estimated values 





(c) L-band. 

Figure 27.- Concluded. 
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